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ABSTRACT 

We present a statistical analysis of the largest X-ray survey of nearby spiral galax- 
ies in which diffuse emission has been separated from discrete source contributions. 
Regression and rank-order correlation analyses are used to compare X-ray properties 
such as total, source and diffuse luminosities, and diffuse emission temperature, with 
a variety of physical and multi-wavelength properties, such as galaxy mass, type and 
activity, and optical and infrared luminosity. 

The results are discussed in terms of the way in which hot gas and discrete X-ray 
sources scale with the mass and activity of galaxies, and with the star formation rate. 
We find that the X-ray properties of starburst galaxies are dependent primarily on 
their star-forming activity, whilst for more quiescent galaxies, galaxy mass is the more 
important parameter. One of the most intriguing results is the tight linear scaling 
between far-infrared and diffuse X-ray luminosity across the sample, even though the 
hot gas changes from a hydrostatic corona to a free wind across the activity range 
sampled here. 

Key words: Surveys - galaxies: spiral - galaxies: starburst - galaxies: ISM - galaxies: 
fundamental parameters - X-rays: galaxies 
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1 INTRODUCTION 

Though it was the Emstem satellite that really began the 
study of the X-ray emission from nearby galaxies, its instru- 
mentation only enabled detailed studies of the very closest 
systems (e.g. M31, M33 and the Magellanic Clouds). It was 
only with the launch of ROSAT (Triimper 1992) that de- 
tailed work on more distant systems, or any reliable spectral 
work, became possible. 

Results from both Emstein (see e.g. the X-ray cat- 
alogue of Fabbiano, Kim & Trinchieri (1992)) and ROSAT 
(e.g. Read, Ponman & Strickland 1997) suggest that normal 
spiral galaxies emit X-rays in the soft X-ray (0.1— 2.4keV) 
band with luminosities of ~ 10^* — 10*^ erg s~^. This emis- 
sion originates from a combination of the integrated output 
of supernova remnants and accreting close binaries (low- and 
high-mass X-ray binaries, cataclysmic variables etc.), and 
hot gaseous diffuse emission. The latter may take the form 
of hot outflows (or winds), as seen emanating from the discs 
of nearby starburst galaxies (Watson, Stanger & Griffiths 
1984, Fabbiano 1988), or hot coronal phases of the ISM, 
seen within more normal galaxies (e.g. Bregman & Pildis 
1994; Wang et al. 1995). So, the X-ray emission from spiral 
galaxies is generally composite and complex. Hence simple 
integrated analyses, which combine the different components 



and extract mean parameters for the whole, are of rather 
limited value. 

In Read, Ponman & Strickland (1997), hereafter Pa- 
per 1, we reported the results of the largest X-ray survey 
of nearby spiral galaxies performed with the ROSAT PSPC 
which involves a careful separation of diffuse and point-like 
emission - a uniform analysis of a sample of 17 nearby spi- 
ral systems covering a wide range in activity and inclination. 
The separation of point sources from the diffuse emission al- 
lows an insight into the properties of X-ray sources within 
external galaxies, and enabled us to establish the prevalence 
and properties of hot galactic winds and halos. We discussed 
the X-ray properties of both the point-source emission and 
the diffuse emission within each of the 17 systems in de- 
tail, together with the general X-ray properties of the whole 
survey. Physical parameters for the hot gas observed within 
many of these systems were also derived. 

We concluded from this study that nearby non-active 
spiral galaxies have (0.1— 2.0keV) X-ray luminosities in the 
range 10'^*-10*^ erg s~^ and that significant amounts (up to 
10^ M0) of diffuse hot (1 - 8 x 10** K) gas are evident in 
most of the systems observed (the amount of gas appear- 
ing to increase with Lx)- This hot gas extends above the 
planes of high-activity galaxies in the form of winds and ha- 
los. Point sources within these galaxies were found to have 
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(0.1 — 2.0 keV) X-ray luminosities in the range of a few x 10''^ 
to a few X 10"^^ erg s~^, and dominate the emission from the 
low-activity systems. Most of the sources with Lx > 10^^ erg 
s~^ were nuclear sources with soft (fcr<,lkeV) spectra. 

In recent years, several related studies have appeared. 
The most similar to our own work is the X-ray mini-survey of 
nearby edge-on spiral galaxies by Dahlem, Weaver & Heck- 
man (1998). This study has several galaxies in common with 
our sample, and there is in general good agreement between 
the two studies for these systems (note that the Dahlem et 
al. (1998) study makes use also of ASCA data and more 
complex spectral fitting). Other recent galaxy surveys using 
ROSAT data include a multivariate statistical analysis of 
integrated spiral galaxy luminosities (Fabbiano & Shapley 
2000), and a ROSAT HRI survey of bright nearby galax- 
ies (Roberts & Warwick 2000) focusing on the point source 
content. 

In the present paper we use the results presented in 
Paper 1 to analyse the statistical properties of the sample 
as a whole. This study is the best of its kind to date, in terms 
of the size of the sample coupled with the detail of the X- 
ray analysis. The way in which the X-ray properties of these 
spiral galaxies, such as their X-ray luminosity, their diffuse 
fraction and the temperature of the diffuse emission, vary 
with other properties, such as galaxy morphology, activity 
and optical and infrared luminosity, are discussed. 

The plan of this paper is as follows. After some brief 
comments (Section Bh on the selection of the original sam- 
ple, the observations and the data reduction methods used, 
the results of statistical properties of the sample as a whole 
are presented (Section 0). Important implications of these 
results are discussed in Section 0, and finally in Section H 
we present our conclusions. 



2 THE SAMPLE, OBSERVATIONS AND DATA 
REDUCTION 

A brief description of the sample and the reasoning behind 
its selection is given below - this is described in greater 
detail in Paper 1. 

The sample comprises 17 nearby 'normal' (i.e. systems 
without strong AGN) spiral galaxies chosen to have one 
(or more) of the following properties: (i) They are X-ray 
bright, having an Einstein (0.2— 4.0keV) X-ray flux (Fab- 
biano et al. 1992) greater than 45 x 10"^^ erg cm~^ s~^, 
(ii) they are infrared bright, having a 100 /xm infrared flux 
(Soifer, Boehmer & Neugebauer 1989) greater than 64 Jy, 
and/or (iii) they have a large value of 'supernova flux' (i.e. 
fsN = LsN/'iT^d^), where Lsn is the estimated supernova 
luminosity and d is the distance) (see Read & Ponman 1995) 
(note that many of these galaxies are also found in the EUV 
sample of nearby spirals of Read & Ponman (1995)). 

Basic properties of the sample galaxies are shown in 
Table qI. For the sake of consistency, the galaxy distances in 
Tully (1988), which are based on Hq = 75km s"^ Mpc"\ 
and assume a Galaxy retardation of 300 km s"'^ from uni- 
versal expansion by the mass of the Virgo cluster, have been 
used. Axial ratios and major axis diameters are taken from 
the Second Reference Catalog of Bright Galaxies (RCBG) 
(de Vaucouleurs, de Vaucouleurs & Corwin 1976). Optical 
(B) luminosities for all the galaxies in the sample are taken 



from Tully (1988) and FIR luminosities are calculated from 
IRAS 60 and 100 ^m fluxes (Soifer et al. 1989; Rice et al. 
1988) using the expression 

Lfir = 3.65 x W [2.58560m™ + Siooi^m] D Lq, 

(e.g. Devereux & Eales 1989), where D is the distance 
(in Mpc) and Seo/j™ a^nd SiooMm , the IRAS 60 and 100 /xm 
fluxes (in Janskys). An asterisk by the galaxy name indi- 
cates a 'starburst' galaxy, i.e. one deflned as here with 
Lfir > O.SSLb- Also given in Table hi is the FIR colour tem- 
perature Smtim/Sioofim (a good measure of the activity), the 
galaxy type (i.e. the mean numerical index of stage along 
the Hubble sequence), from the RCBG, and, where applica- 
ble, the Einstein (0.2—4.0 keV) X-ray luminosity, calculated, 
for the distances tabulated, using the fluxes given in Fab- 
biano et al. (1992), Trinchieri et al. (1988) and Trinchieri 
et al. (1990). For convenience, the ROSATX-iay luminosity 
(as given in Table H) is also given. The values of the Galactic 
hydrogen column are taken from the neutral hydrogen radio 
results of Stark et al. (1992). 

Observations, data reduction and results for the indi- 
vidual systems are described in full in Paper 1. However, a 
brief summary of the reduction techniques is given here for 
convenience. 

AU PSPC datasets (obtained from the UK ROSATData 
Archive Centre at the Department of Physics and Astron- 
omy, Leicester University, UK) were reduced in exactly the 
same way using the Starlink ASTERIX X-ray analysis sys- 
tem. 

In each case, the data were 'cleaned' of high-background 
times, and the point source search program PSS (Allan, 
Ponman & Jeffries, ASTERIX User Note 4) was used to 
search for point-like sources, the positions of which were then 
cross-correlated with a variety of stellar and non-stellar cat- 
alogues. An integrated (source -I- unresolved emission) spec- 
trum was formed from a circular area of diameter slightly 
greater than the optical diameter of each galaxy (so as to in- 
clude any diffuse emission above the plane in edge-on galax- 
ies). 

A key aspect of the analysis performed in Paper 1, was 
the separation of the emission into diffuse and source compo- 
nents. In order to separate the diffuse from the source emis- 
sion, data were removed at the position of each PSS source, 
and the remaining data were then collapsed into a spectrum 
and corrected for vignetting effects and exposure time. To 
account for diffuse flux lost in the source- removal procedure, 
the diffuse spectrum was renormalised using a 'patched' im- 
age, where the holes left after source removal were flUed 
by bilinear interpolation. In Paper 1 we discussed tests of 
the reliability of these procedures for separating pointlike 
and diffuse emission. As might be expected, we found a de- 
pendence on the angular size of galaxies, and also on their 
activity. The worst case is distant galaxies with strong dif- 
fuse emission. In this case a signiflcant fraction of the diffuse 
flux may be erroneously attributed to sources. In the present 
sample, we estimate that at most three systems could be 
affected at a level as large as 30% of their diffuse flux. It 
should also be stressed that we have covered a large range 
in Hubble galaxy type, from galaxies with large elliptical- 
like bulges (Sa spirals) to galaxies with very small bulges 
and dominating spiral arms (Sd spirals). Furthermore, al- 
though the survey galaxies have been chosen as containing 
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Table 1. Physical and multiwavelength properties of the survey sample (see text for details). 



no dominating AGN, a modest, low-level range in AGN ac- 
tivity does exist across the survey. 

Source spectra were binned directly from the raw data, 
though in order to arrive at a true source spectrum how- 
ever, it was necessary to remove the flux from both the local 
true background and the local diffuse emission, and this re- 
quired quite an intricate procedure (described in Paper 1). 
All source spectra, integrated spectra and diffuse emission 
spectra were fitted with standard spectral models. Power 
law and thermal bremsstrahlung models were fitted to each 
of the point source spectra, whilst power law, Raymond & 
Smith (1977) hot plasma, and differential emission measure 
models (Raymond and Smith models integrated over a range 
of temperatures) were fitted to the diffuse and integrated 
spectra. Minimum chi-squared fitting was used for the dif- 
fuse and integrated spectra, and so a goodness of fit could 
be obtained. For the point source spectra, which generally 
contained far fewer counts, Gaussian statistics cannot be 
assumed, so a maximum likelihood criterion was used. 

Assuming the diffuse emission to be hot gas, basic phys- 
ical properties of this gas could be inferred on the basis of as- 
sumptions about the geometry of the emission. In the bubble 
model, where the gas is assumed to be contained in a spher- 
ical bubble of radius r, the fitted emission measure rjn^V 
(where rj is the 'filling factor' - the fraction of the total vol- 
ume V which is occupied by the emitting gas) can be used 
to infer the mean electron density Tie, and hence the total 
mass Mgas (oc neV), thermal energy E^.-^ (oc n^VT, where 
T is the temperature), and cooling time t^oQi (c>c E^^^^L^^) 
of the hot gas. 

Table presents a summary of the results obtained in 
Paper 1. After the galaxy name (asterisks again indicating 
the 'starburst' galaxies), is given the integrated (i.e. source 
plus diffuse) (0.1— 2.0keV) X-ray luminosity (also given in 
Table pi). The next two columns give the number of sources 
detected within the galaxy and the integrated source X-ray 



luminosity. The remaining eight columns deal with the dif- 
fuse emission. We list the X-ray luminosity of the diffuse 
emission, the 'diffuse fraction' (the ratio of Lx to Lx 
plus I/x*^), and the fitted temperature of the diffuse emis- 
sion. The last five columns deal with the 'bubble' model gas 
parameters: the radius r of the assumed bubble is taken as 
the average radius of the lowest contour level seen in each 
diffuse emission image in Paper 1; this is followed by the 
mean electron density, total mass, thermal energy and cool- 
ing time of the hot gas. 



3 RESULTS 

Correlation tests were performed between almost all of pos- 
sible pairs of the parameters listed in tables hj and Q, and for 
a variety of galaxy selections. The most interesting results 
of regression and Spearman rank correlation analyses per- 
formed on the data are given in Table H as follows: The Y 
and X parameters (cols. 1 & 2), the galaxy selection criteria 
used (whether all the galaxies, just the starbursts [S], just 
the normal galaxies [N], or any individual omissions [e.g. 
no M82]) (col. 3), and the values and associated Icr errors 
for the regression gradient m (col. 4) and constant c (col. 5). 
Here the assumed regression is of the form Y = mX + c. 
Also given (col. 6) , is the Spearman rank-order correlation 
coefficient rg (lying between -1 and -1-1), and (col. 7) the 
associated significance of rs, Tg. Here, Ts is distributed ap- 
proximately as Student's distribution with N — 2 degrees of 
freedom. Only fit results for which the regression software 
produced a stable solution are listed in Table H. For exam- 
ple, for the log Lx '■ log Lb relation, no fit for the starburst 
subsample could be derived, though a stable solution was 
obtained when M82 was omitted from the sample, and this 
is listed. Further results, for instance Spearman rank correla- 
tion coefficients for pairs of parameters where no regression 
fits were found, are given in the following text. 
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Table 2. Physical properties derived from the X-ray analysis, summarised from Paper 1, including integrated (total) and point source 
and diffuse luminosities, number of X-ray sources detected, and a variety of inferred properties of the hot gas (see text for details). 



From amongst the large number of available relations, 
we concentrate below on those which are available here for 
the first time, as a result of our analysis: the individual 
source and diffuse contributions to the X-ray emission, and 
properties derived from these. We examine the variation of 
these components with galaxy mass and activity, discuss the 
relation to emission in other wavebands, and finally investi- 
gate some of the derived gas parameters, and the variation 
in X-ray properties with Hubble galaxy type. 

Scatter plots follow a standard format: individual 
galaxy names are given, centred on their particular positions 
in the relevant parameter space (brackets indicate that the 
parameter position of the galaxy lies closer to the nearest 
non- bracketed galaxy than is reasonable to plot). Starburst 
galaxies are shown in a bolder font and underlined. Where 
relevant, the regression fits from Table 0, for the full sample 
(bold line), and for the (S)tarburst and (N)ormal subsam- 
ples (dashed lines) are also plotted. 

3.1 Mass and activity 

Before proceeding further, we need to clarify how we define 
the 'mass' and 'activity' of our galaxies. We take the optical 
luminosity. Lb, as our measure of mass. This is not ideal, 
since the blue luminosity to mass ratio is quite sensitive 
to the age of a stellar population. Luminosity in the near 
infrared would be a more robust indicator of stellar mass, 
but is not available for many of the galaxies in our sample. 
Several studies over the years (e.g. Devereux & Young 1991; 
Joseph et al. 1984) have linked the far-infrared luminosity 
of a galaxy to the star- formation rate (SFR), and hence 
(for a given initial mass function) to the supernova rate. 
We therefore take these rates to be proportional to Lfir- 
Therefore, we can also consider the 'activity', i.e. the 'SFR 
per unit mass', to be measured by Lfir/ Lb- 



Given the parameters at our disposal, we could consider 
other indicators of mass and activity. For example, the far- 
infrared colour temperature SeoMm/S'iooum (see e.g. Telesco 
et al. 1988) is an attractive activity indicator, except that 
it has rather large statistical errors for many of our galax- 
ies. An alternative possible measure for the mass might be 
galaxy volume, defined as (diameter'^ x axis-ratio) ^^ Com- 
paring these alternative measures, it is encouraging to note 
that they are well correlated, as can be seen in Fig. 1 (top 
and middle) - Lb is strongly correlated to galaxy volume 
{Ts — 3.2), and Lfir/ Lb with the far- infrared colour tem- 
perature (Ts = 2.5). This strengthens our confidence that 
the chosen parameters are good tracers of mass and activity. 

Since we wish to separate the effects of mass and activ- 
ity on galaxy properties, it will help a great deal if mass and 
activity are independent variables for our galaxy sample. If, 
for example, large galaxies in the sample tended to have 
higher activity than small ones, it would be difficult to dis- 
entangle the effects of the two variables. Fortunately, there 
is no strong correlation between the mass tracer and the 
activity tracer (see Fig. 1 (bottom)). Our sample does ap- 
pear to be lacking in high-LFiR-iow-Ls systems, i.e. high- 
activity dwarfs. This is very likely due to them being rare, 
and hence distant, and thus not meeting the flux limits set 
for the sample. Despite this, the Spearman rank significance, 
Ts, for any correlation between our mass and activity pa- 
rameters lies well below 1 not only for the total sample, but 
also for the normal and starburst subsamples. This makes 
the analysis of the following results a good deal easier. 



3.2 Variations in X-ray properties with mass and 
activity 

First, we compared the different components of the X-ray 
emission with the galaxy mass for the total sample, and 
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-24.87±0.38 


0.60 


2.9 






S, no M82 


[7] 


1.11±0.39 


-8.49±0.24 


0.70 


2.2 






N[9] 




1.44±0.12 


-23.26±0.19 


0.90 


5.6 


logLdiff 


log Lb 


All [17] 




1.89±0.25 


-43.06±0.48 


0.49 


2.2 






S, no M82 


[7] 


0.92±0.42 


-0.33±0.26 


0.66 


2.0 






N[9] 




1.91±0.14 


-43.88±0.22 


0.95 


8.2 


logL^Jc 


log Lb 


All [17] 




1.23±0.18 


-14.05±0.34 


0.64 


3.2 






S, no M82 


[7] 


1.37±0.43 


-20.34±0.26 


0.63 


1.8 






N[9] 




1.15±0.13 


-10.72±0.19 


0.76 


3.1 


log Lpjji 


log Lb 


All [17] 




1.48±0.30 


-21.63±0.56 


0.51 


2.3 






S, no M82 


[7] 


0.92±0.29 


3.13±0.18 


0.80 


3.0 






N[9] 




1.22±0.28 


-10.63±0.42 


0.71 


2.7 


logLdiff/ifl 


log Lb 


All [17] 




0.89±0.25 


-43.09±0.48 


0.50 


2.2 






N[9] 




0.91±0.14 


-43.88±0.22 


0.87 


4.7 


logLdiff/La 


log Lfir/ Lb 


All [17] 




0.75±0.19 


-3.70±0.45 


0.71 


3.9 






S[8] 




0.96±0.23 


-3.59±0.23 


0.58 


1.8 


log L^^y Lb 


logLpjR/Ls 


All [17] 




0.22±0.05 


-4.28±0.23 


0.55 


2.6 


logr-difj 


log Lb 


[14] 




0.65±0.12 


-27.51±0.16 


0.66 


3.0 






S[7] 




0.28±0.10 


-11.08±0.07 


0.80 


3.0 






N[7] 




0.76±0.20 


-32.34±0.20 


0.63 


1.8 


/diff 


log Lb 


All [17] 




0.32±0.09 


-13.41±0.16 


0.43 


1.8 






N[9] 




0.36±0.09 


-15.14±0.13 


0.82 


3.8 


logidiff 


logoff 


[14] 




2.5±1.6 


40.84±0.69 


0.48 


1.9 


log Mgas 


log LpiR 


[14] 




0.77±0.16 


-25.24±0.47 


0.70 


3.4 






S[7] 




1.05±0.43 


-37.34±0.30 


0.70 


2.2 






N[7] 




0.94±0.35 


-32.44±0.63 


0.40 


1.0 


log Lx/ Lb 


Type 


no M82 [16] 


-0.075±0.035 


-3.48±0.31 


-0.44 


-1.8 






N[9] 




-0.083±0.028 


-3.63±0.21 


-0.71 


-2.7 



Table 3. Main results of regression analyses (assuming slopes of the form Y = mX + c, unweighted on Y and X) performed on the 
galaxy sample. The Table lists the Y and X parameters, the galaxy selection criteria used, and the values and associated la errors for the 
regression gradient m and the intercept c. Also given is the Spearman rank-order correlation coefficient rg and its associated significance 
Tg (see text for details). 



for the normal and starburst subsamples. Fig. 2 shows the 
variation in the difTuse X-ray emission fraction (/dig) with 
the 'mass' (Lb)- As can be seen in the figure (and in Table 
bl), a gradual increase in the diffuse X-ray emission fraction 
with mass is observed, i.e. high mass galaxies have a smaller 
fraction of their X-ray emission in discrete sources. This re- 
lationship is driven by the normal subsample, and is not seen 
in the starbursts. 

Comparing source and diffuse X-ray luminosity per unit 
galaxy mass (i.e. L^^/Lb and Lx /Lb) with galaxy mass, 
we find no significant trend in source luminosity (per unit 
galaxy mass). However, the production of diffuse X-ray emis- 
sion per unit galaxy mass is seen to correlate strongly with 
galaxy mass (see Table H) for the normal galaxy subsample 
(Ts = 4.5). 

The same X-ray parameters were compared with the ac- 



tivity tracer, Lfir/ Lb- Diffuse fraction (fdis) is seen to rise 
with activity within the starburst sample {Ts — 2.4), but 
shows no significant trend within the normal galaxy sample. 
Diffuse X-ray luminosity per unit galaxy mass {Lx /Lb), 
Fig. 3 (top), shows a trend with activity across the whole 
galaxy sample. Starbursts have (by definition) high activity, 
and also higher diffuse X-ray emission than normal galaxies. 

The source X-ray luminosity per unit galaxy mass (see 
Fig. 3 (bottom)) does show a general rise with activity from 
normal to starburst galaxies within the total sample {Ts — 
2.6), though the trend is not significant within the individual 
subsamples {Ts < 1 in both cases). 

The broad picture which emerges from the above is that 
normal galaxies appear to be very much more influenced by 
their mass than by their activity. For starburst galaxies, on 
the other hand, it is activity that is the dominating factor. 
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Figure 1. 

Relationships between (top) the two suggested 'mass' parameters, 
Lg and 'volume', (middle) the two suggested 'activity' param- 
eters, Lpir/ Lb and the far-infrared colour temperature, and 
(bottom) the adopted mass and activity parameters, Lg and 
Lfir/Lb- 



Figure 2. 

The variation in the diffuse X-ray emission fraction (/(jjff) with 
galaxy 'mass' {Lb)- 



We will discuss the properties of the diffuse X-ray emis- 
sion in detail later, but comment here on some of the more 
noteworthy effects of galaxy mass and activity on the proper- 
ties of the gas. The derived temperatures of the diffuse X-ray 
emission shows no really significant trends with either mass 
or activity. Only a low-significance positive correlation in 
diffuse emission temperature with activity is observed, and 
then only in the starburst subsample {Ts = 1.3). The radial 
extent of the diffuse emission does show significant positive 
correlations with galaxy mass. Fig. 4 shows the relationship 
between the log of the radial extent of the diffuse emission 
(r; see Table 0) and the galaxy mass tracer Lb- As can be 
seen from the figure, and from Table M, a positive trend is 
seen, both in the total sample, and in the individual star- 
burst and normal subsamples. Larger galaxies posess larger 
diffuse emission features. 

Surprisingly, only a low- significance positive correlation 
is seen between the radial extent of the diffuse emission and 
galaxy activity. For the total sample, only a Spearman rank 
coefficient of 0.42 (corresponding to Ts: 1.6) is obtained. 
Nothing significant is seen within the individual normal and 
starburst galaxy subsamples. 



3.3 X-ray— Far-Infrared luminosity relationships 

Since Lfir is a measure of star formation rate, a posi- 
tive relation with X-ray emission is to be expected. Fig. 5 
shows the relationships between the various components of 
the X-ray emission and the far-infrared luminosity. ROSAT 
(0.1— 2.0keV) luminosities of the total (top), the diffuse 
(middle) and the source (bottom) emission components are 
shown plotted against the far-infrared luminosity for the sur- 
vey sample. 

Lx correlates tightly with Lfir - the scatter of the 
galaxies about this relationship is remarkably small, the two 
main outliers being NGC 4258 and NGC 1291. That these 
two galaxies appear different from the norm is very likely due 
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Figure 3. 



The variation of (top) the diffuse X-ray luminosity per unit galaxy 

/Lb) and (bottom) the source X-ray luminosity per 
rsrc 



mass (L'^^l 

unit galaxy mass (L^'^ /Lg) with the activity (Lpir/Lb) for the 



galaxy sample. 

to the facts that (i) NGC 4258 is the most AGN-dominated 
system in the survey, and is known to have excess diffuse 
X-ray emission associated with two central AGN-produced 
'jets' (Paper 1; Pietsch et al. 1994), and (ii) NGC 1291 is of 
type Sa, and is the earliest type system in the survey, with a 
very elliptical-like bulge and a consequently low far-infrared 
luminosity (Bregman, Hogg & Roberts 1995). The slope of 
the Fig. 5 (top) relationship {Lx oc I/j^f^j) appears to be 
less than unity, and no significant change in slope is seen 
between the normal and starburst galaxies (see Table H) . 

The relationship between the purely diffuse X-ray lu- 
minosity and the far-infrared luminosity (Fig. 5 [middle]) 
appears very similar to that of Lx-Lfir- Lx increases 
with Lfih, and very little scatter is seen, apart again from 
NGC 4258 and NGC 1291. However, here the slope is unity, 
Lx oc Lfir- Again, no significant change is seen between 
the normal and the starburst galaxies. The corresponding 



Figure 4. 

The relationship between the log of the radial extent of the diffuse 
emission (r; see Table bl) and the galaxy mass tracer Lg for the 
total sample. 



correlation coefficients are generally even higher than for 
Lx-Lfir^ especially for the starbursts. 

If the total X-ray luminosity is increasing less rapidly 
than the far-infrared luminosity, while the diffuse X-ray com- 
ponent keeps pace, then this suggests that it must be the 
remaining component, the source X-ray component, is lag- 
ging Lfir- Indeed, as is seen in Fig. 5, (bottom), and in 
table M, the Lx{src)-LFiR slope is flatter (Lx(src)oc L^f^). 
Note again, the lack of any significant difference between the 
normals and the starbursts. 



3.4 X-ray Optical luminosity relationships 

Fig. 6 shows the relationships between the total (top), the 
diffuse (middle) , and the source (bottom) X-ray luminosity, 
with the optical luminosity (a measure, remember, of the 
mass of the galaxy) . 

As expected, Lx, Lx and L^'^ all increase with Lb- 
The total X-ray:optical slope is consistent with the recent 
Emstein work of Fabbiano & Shapley (2000). Little scatter 
is seen at the low Lb end of each plot, whereas a good deal of 
scatter is seen at higher Lb - This is evident in the correlation 
coefficients and regression fit results in Table 0. In both the 
total and diffuse cases, the normal galaxy subsamples are 
very well correlated, whereas stable fits (and even then, with 
a great deal of scatter) are only obtained for the starburst 
subsample when M82 is removed. Note that while for the 
starbursts, Lx appears well correlated with Lfir, for the 
normals it is very well correlated with Lb- 

In the total and diffuse cases, several things can be seen. 
Firstly, for the whole sample, Lx (and Lx ) rise faster than 
Lb, with slopes significantly greater than unity. Secondly, 
in both cases, the total sample slopes are essentially identi- 
cal to their respective normal galaxy subset slopes. Thirdly, 
though no regression fits are obtainable (in either case) for 
the full starburst galaxy subsets, when the 'irregular' M82 
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Figure 5. 

ROSAT (0.1-2.0 keV) total (top), diffuse (middle) and source 
(bottom) X-ray luminosity versus far-infrared luminosity for the 
survey sample. Starburst galaxies are marked in a bold font and 
underlined. Regression fit lines to the total sample (bold line) 
and to the (S)tarburst and (N)ormal subsamples (dashed line) 



Figure 6. 

ROSAT (0.1-2.0 kcV) total (top), diffuse (middle) and source 
(bottom) X-ray luminosity versus optical luminosity for the sur- 
vey sample. Starburst galaxies are marked in a bold font and 
underlined. Regression fit lines to the total sample (bold line) 
and to the (S)tarburst and (N)ormal subsamples (dashed line) 
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is removed from the analysis, the regression slopes obtained 
are significantly lower, more like unity in both cases, and are 
significantly difi'erent from both the total and normal subset 
slopes. 

For the source emission, the slope of L^x^:Lb is lower 
than for Lx '-Lb for the total sample, and the starburst and 
normal subsamples - all being nearly consistent with unity. 
In contrast to the situation with diffuse emission, there is 
no indication of any flattening of the relation for starbursts. 



3.5 The diffuse emission and hot gas properties 

Fig. 7 shows the relationship between the X-ray luminos- 
ity, Lx , and temperature, T^iQ, of the diffuse emission. 
The L : T relation has been a topic of great interest in 
the study of X-ray emission from galaxy clusters, since it 
is easy to show that for hot gas with universal mean den- 
sity, filling potentials which scale in a self-similar way, Lx 
is expected to scale as T^'^A(T), where A(T) is the tem- 
perature dependence of the X-ray emissivity. In the case of 
bremsstrahlung, A(r) oc T^", so one expects Lx oc T^. For 
gas with T < 2keV, such as we observe in spiral galaxies, 
A(r) is very flat, and therefore L would rise more slowly 
that T^. 

Whilst self-similar potentials and universal gas densities 
may be reasonable first order expectations in the case of 
clusters, this is much less obviously the case for individual 
spiral galaxies. Here the additional physics of heating and 
cooling which affect the baryons, may be expected to modify 
both the potentials, and especially the gas fractions and gas 
density profiles. In addition, we know from its morphology 
(Paper 1) that the diffuse X-ray emission in some systems 
is dominated by a starburst-driven wind, whilst in others 
it arises from a gravitationally bound corona, or from hot 
bubbles within the disc. 

We should not be surprised, therefore, to see a large de- 
gree of scatter in the relation shown in Fig. 7. Nonetheless 
there is a marginally significant trend {Ts — 1.9) with a log- 
arithmic slope of 2.5 ± 1.6. The correlation does not result 
from a scaling with galaxy size, since T^^ff does not corre- 
late with Lb ■ Rather the relation appears to arise principally 
from the fact that non-starbursts tend to have both lower 
temperatures and lower diffuse X-ray luminosities than star- 
burst galaxies. 

The relationship between the temperature and the in- 
ferred mass of the hot gas is shown in Fig. 8. There is a 
tendency for more hot gas to be found in those systems 
exhibiting higher gas temperatures. No acceptable regres- 
sion fits could be obtained, but the correlation coefficients 
[total sample: 0.57 {Ts 2.4) starbursts: 0.87 {Ts 3.9) nor- 
mals 0.77 {Ts 1.7)] conflrm what can be seen directly from 
Fig. 8 - the correlation arises principally from the starburst 
galaxies. In a similar vein Fig. 9 shows that there is a rela- 
tionship between the inferred mass of the hot gas with the 
far-infrared luminosity, which is especially marked amongst 
the starbursts. 
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Figure 7. 
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Log-log plot of the diffuse emission temperature TJ^jff against the 
diffuse X-ray luminosity (0.1— 2.0keV). 
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Figure 8. 

The relationship between the diffuse gas mass and the diffuse gas 
X-ray temperature for those members of the sample for whieh 
reliable measurements have been obtained. 



3.6 Variations in X-ray properties with Hubble 
galaxy type 

Fig. 10 shows the relationship between the diffuse gas tem- 
perature and the Hubble type of the galaxy. Although we are 
dealing with small numbers here, there is a suggestion, for 
the normal galaxies, that the earlier Hubble types contain 
hotter gas (correlation coefficient -0.54, Ts -1-4), whereas 
for the starburst galaxies, the later type systems are hotter 
(correlation coefficient 0.68, Ts 2.1). 

The relationship between X-ray to optical luminosity 
ratio and Hubble type is given in Fig. 11. Again for normal 
galaxies, a weak negative correlation is seen (i.e. as we move 
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Figure 9. 

The log-log relationship between the diffuse gas mass and the far- 
infrared luminosity for those members of the sample for which 
reliable measurements have been obtained. 



Figure 11. 

The relationship between the X-ray to optical luminosity ratio 
and the galaxy Hubble galaxy type for sample members (M82 
has been marked with asterisks, as it is an irregular galaxy). 
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Figure 10. 

The relationship between the diffuse gas temperature and the 
galaxy Hubble galaxy type for those members of the sample for 
which reliable measurements have been obtained (M82 has been 
marked with asterisks, as it is an irregular galaxy). 

towards earlier Hubble type, the X-ray to optical luminosity 
ratio appears to increase). However, nothing concrete can 
be said here regarding the starburst sample. 

Finally, it is worth noting the case of the earliest type 
system in the survey, NGC 1291. This Sa-type galaxy ap- 
pears to have a normal X-ray to optical luminosity ratio 
(Fig. 11, Fig. 6), but a very high X-ray to far- infrared lu- 
minosity ratio (Fig. 5). It is known that this system's X- 
ray luminosity and radial surface brightness distribution are 
like that of elliptical galaxies, and that, like ellipticals, the 
log Lfir — log Lb ratio is small (Bregman et al. 1995). This 



can be seen in Figs. 1 & 3. Throughout this paper, NGC 1291 
is seen to appear elliptical-like in that its far-infrared lumi- 
nosity is reduced in comparison with its X-ray and optical 
luminosity, whereas the next earliest galaxies, NGC 891 and 
NGC 3628, appear like typical spirals. 



4 DISCUSSION 

In the previous section, we presented the main results from a 
statistical analysis of the galaxy sample. In discussing these 
results, we concentrate here on the parameters which have 
not generally been available from earlier studies. In partic- 
ular, we have been able to split each system's X-ray emis- 
sion into a source component and a diffuse, gaseous compo- 
nent, and also to investigate the spectral properties of the 
hot gas. Our 'new' X-ray parameters are therefore the total, 
source and diffuse X-ray luminosities, and the diffuse emis- 
sion temperatures (together with related parameters, such 
as gas mass etc.). 

4.1 X-ray emission as a function of mass and 
activity 

The bro ad p icture which emerges from the results presented 
in sect. 



3.2 



is that for starburst galaxies. X-ray properties 
scale with activity rather than mass. For the normal galax- 
ies, in contrast. X-ray properties are strongly correlated with 
mass, and not with activity. 

Large normal galaxies, irrespective of activity, appear 
to possess more diffuse emission (here assumed to be pri- 
marily hot gas - see Section 2 and Paper 1), than smaller 
galaxies. This is not a simple a matter of scaling with size, 
since the diffuse X-ray luminosity per unit galaxy mass, and 
the diffuse emission fraction both increase with Lb- Since 
the luminosity of hot gas scales as n^V, a simple size scal- 
ing leads to Lx oc Lb, if mean density and gas fraction 
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do not change with galaxy mass. It appears therefore that 
large galaxies are able to heat or retain a larger fraction of 
their gas than do smaller galaxies, or they compress it to 
higher mean densities. One plausible picture is that hot gas 
from active star formation regions is able to escape more 
readily from the shallower potential wells of small galaxies, 
systematically lowering their hot gas fraction. 

In the case of starburst galaxies, the fraction of diffuse 
emission (i.e. hot gas) is primarily determined by activity 
rather than galaxy mass. This suggests that the influence 
of the galaxy potential, important for the normal galaxies 
above, is no longer as important. The activity now domi- 
nates, having reached such a level that whether the galactic 
potential is large or small has little effect. 

This is all consistent with the idea, supported by the 
morphology of the emission (Paper 1), and by simulations 
and observations at other wavelengths, that in normal galax- 
ies one sees bound 'coronal' gas, whilst in starburst galaxies 
the bulk of the hot gas is a freely escaping hot wind. 

Turning to the source component of the X-ray emission, 
this is found to scale approximately linearly with Lb for 
both normal and starburst galaxies. This result is readily 
understood if the discrete source complement scales with 
the overall stellar mass, and suggests that the separation 
of source and diffuse contributions in our sample has been 
successfully achieved. Although the principal scaling is with 
galaxy mass, there is also a trend of increasing L^x'^ /Lb 
with galaxy activity (Fig. 3). Since we believe that source 
confusion does not have a major effect on our separation of 
diffuse and source emission (section 2), the implication is 
that there is a significant rise in the incidence of discrete X- 
ray sources in systems which are actively forming stars - i.e. 
these galaxies contain a significant contribution from short- 
lived sources, such as high mass X-ray binaries [HMXB] or 
young SNR. In contrast, we know that the source population 
in the Milky Way is dominated by the older low mass X-ray 
binary (LMXB) population. 



transition from an essentially hydrostatic hot corona to a 
wind occurs within this range. Furthermore, Read & Pon- 
man (1998) found that, in even higher activity interacting 
and merging systems, a change in the X-ray:FIR slope is 
seen; the logarithmic X-ray:FIR slope is seen to drop signifi- 
cantly to around 0.3—0.4. It was suggested that this relative 
deficit in X-ray emission might be due to much of the in- 
put energy from supernovae and stellar winds in these very 
high-activity systems, being converted into the kinetic en- 
ergy of the huge plumes of ejected gas seen in these systems, 
rather than contributing to the X-ray emission. In contrast, 
the high-activity systems in the present sample do not be- 
have like these merging systems, even though in several cases 
(M82, NGC 253, NGC 3079 etc.), large plumes of hot gas 
do exist. 

The work of Fabbiano & Shapley (2000), who anal- 
yse the results of a multiwavelength statistical study of the 
Emstem sample of spiral and irregular galaxies, is consis- 
tent with the above results. These authors derive X-ray:FIR 
slopes very similar to our present study. They also see a sug- 
gestion of the X-ray:FIR slope flattening, though in their 
data it is still consistent with a constant slope. This is of in- 
terest as the Fabbiano & Shapley (2000) study did include 
very high-activity galaxies, galaxies where, as in the Read 
& Ponman (1998) study, kinetic energy losses had become 
important. 

It is well established (e.g. Read & Ponman 1995) that 
the X-ray luminosity of late-type galaxies is only a small 
fraction of their total supernova luminosity in all cases. The 
bulk of the supernova energy is apparently radiated in the 
FIR. It appears, however, that an approximately constant 
fraction of the energy released during the life and death of 
stars is radiated by hot gas in both normal and in moderate 
starburst galaxies, but the massive kinetic energy losses suf- 
fered in the ultraluminous outbursts associated with galaxy 
merging, reduces the fraction of the energy available for ra- 
diating in the X-ray band. 



4.2 The effect of star formation rate 



We saw in Sect. 3.S that the diffuse X-ray emission rises 
proportionally with the far-infrared luminosity (i.e. SFR), 
whilst source X-ray luminosity rises less steeply. Combin- 
ing diffuse and source contributions gives the observed 
Lx(total)(x L^p^jB relationship. There is therefore a shift in 
the balance between source and diffuse X-ray luminosity as 
the SFR increases. This is quite easy to explain in terms 
of more of the high-SFR sources being diffuse in nature - 
supernova remnants and collections thereof, superbubbles, 
supernova-produced fountains, halos and coronae of hot gas, 
and galactic winds. Furthermore, we would not expect the 
source X-ray luminosity to keep pace with Lfir (i.e. the 
current SFR), as it is dominated by an older component 
(mostly LMXBs). The above explanation is strengthened by 
the fact that the diffuse X-ray component is seen to be more 
tightly correlated with Lfir than the source component . 

We also found that there appears to be no significant 
change in the X-ray:FIR slope (in any component, whether 
total, diffuse or source) as we go from (low-activity) normal 
galaxies to (high-activity) starburst galaxies. The lack of 
any change in the response of the diffuse luminosity to SFR 
is interesting and surprising, given that we believe that the 



4.3 Diffuse components and their spectral 
properties 

Before discussing the spectral results relating to the diffuse 
components, we should consider first what these diffuse com- 
ponents might be. As discussed in Paper 1 and again in Sec- 
tion 2, we believe the diffuse components observed to be due 
mainly to hot gas produced within these systems. It is be- 
lieved, both thro ugh the literature and through the results 



discussed in Sect. 4.1, that two forms of hot X-ray emitting 
gas can exist in the halos of these spiral galaxies. 

Some edge-on systems within our sample (e.g. NGC 891 
and NGC 4631; see Paper 1) exhibit what is believed to 
be coronal gas (Bregman & Pildis 1994; Wang et al. 1995). 
This cool (1—3 million K) gas, is thought to bubble out of 
the plane of the galaxy through galactic 'chimneys' and 
'fountains' (formed by localized high-activity star-forming 
regions within the disc), then to fall ballistically back to 
the plane (Norman & Ikeuchi 1989). This X-ray emitting 
gas is thought to be in approximate hydrostatic equilibrium 
and occupies almost the entire volume. Its filling factor r] 
is therefore thought to be close to unity and this has been 
backed up with numerical multi-gas-phase simulations (e.g. 
Rosen & Bregman 1995). 
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Other edge-on systems within our sample (e.g. 
NGC 253, M82 and NGC 3079; see Paper 1) exhibit what 
is believed to be an outflowing wind. Many of these struc- 
tures appear to be very extended, having temperatures of 
T~4-8xl0'' K. The structure of the gas is not weU 
known. On the one hand, if the gas filling factor rj is close 
to unity, then these features could be due to plumes of gas 
flowing out from the galaxy. On the other hand, simulations 
(Suchkov et al. 1994, Strickland & Stevens 2000) and obser- 
vations (Heckman et al. 1993, Strickland et al. 2000) suggest 
that in these systems, most of the X-ray emission is due to 
small 'clouds', shock-heated by a fast (~ 3000km s~^), hot 
(3> 10^ K) wind, produced through the efficient thermaliza- 
tion of massive stellar winds and supernovae ejecta within 
a starburst nucleus. This wind is believed to be orders of 
magnitude less luminous in the ROSAT band (Suchkov et 
al. 1994) . In this scenario, the filling factor rj of the ROSAT 
X-ray emitting gas is very much smaller than unity. 

Note that if this is the case, the gas masses inferred 
for the wind systems under the assumption r) = 1 may be 
seriously overestimated. The trends in Mgas from quiescent 
to more active galaxies, plotted in Figs. 8 & 9, may therefore 
be misleading. It is likely that the gas mass estimates for 
coronal systems are reasonably representative, whilst those 
for winds may be too high. 

Considering edge-on systems, in which the morphol- 
ogy of the hot gas can be most clearly seen, we note that 
the diffuse emission in the 'coronal' systems (NGC 4631 
and NGC 891) appear cooler than those seen within the 
known 'wind' systems (NGC 253, M82 and NGC 3079). As 
suggested in Paper 1, this temperature diagnostic suggests 
that the diffuse emission structures seen in NGC 5005 and 
NGC 5055 are more likely to be due to winds than coronae. 

Diffuse structures are also seen in many of the face-on 
systems (e.g. M33, M51, MlOl). We have no idea on ge- 
ometrical grounds, how far these diffuse emission features 
extend above the plane. As discussed in Paper 1, the dif- 
fuse emission seen in M33 is low-power, and only occupies 
the central area of the disc. In M51 and MlOl however, the 
diffuse features observed are both large, powerful, and fairly 
uniform. Their respective positions in Fig. 8 suggest, assum- 
ing the coronal/wind argument described above, that the 
diffuse emission in M51 is dominated by galactic wind emis- 
sion, whereas in MlOl, it is dominated by coronal emission. 

Presumably, there is nothing to stop both coronae 
and winds existing simultaneously in a given galaxy, and 
this may have been observed in the merging galaxy pair 
NGC 4038/9 (Read, Ponman & Wolstencroft 1995). 



5 CONCLUSIONS 

We have reported the results of a statistical analysis of a 
study of 17 spiral galaxies with the ROSAT PSPC. Con- 
tributions to the X-ray emission from discrete sources and 
diffuse emission have been carefully separated, allowing us 
to probe both the hot gas and X-ray sources. The most in- 
teresting conclusions from this work are: 

• In general, the X-ray properties of normal galaxies are 
governed by the galaxy size, whereas the X-ray properties 
of starburst galaxies are governed by the galaxy activity. 



• Larger normal galaxies contain more hot gas per unit 
mass than smaller ones. We interpret this as resulting from 
better retention of hot gas in the deeper potential wells of 
larger galaxies. 

• In contrast, the diffuse emission per unit mass from 
starburst galaxies scales with activity rather than galaxy 
size. This is consistent with the idea that the hot gas in these 
systems is substantially unbound, and its density is therefore 
determined by the rate of hot gas production rather than by 
the retaining potential. 

• The X-ray emission per unit mass from discrete sources 
does not depend on the total galaxy mass, but does rise 
somewhat with galaxy activity. This implies that in the more 
active systems, young objects such as HMXB or compact 
SNR make a significant contribution to the X-ray source 
population. 

• Fits to Lx'-Lb and Lx'-Lfir trends agree well with pre- 
vious studies. Extracting the behaviour of the hot gas, we 
find a remarkably tight linear scaling with Lfir which ap- 
plies across the activity range covered by this study, though 
it appears to be broken in ultraluminous mergers. The rea- 
son for the continuity in this relation across the transition 
from coronal to wind emission is far from clear, and the au- 
thors hope that future hydrodynamical modelling work will 
address this issue. 

• Gas in the 'coronal' systems appears to be cooler than 
that seen in the 'wind' systems. This may help to diagnose 
the state of hot gas in face-on systems. 
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